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I. INTRODUCTION 



Links between two well known methods: methods of zero-range and non-overlapped (muffin- tin) potentials 
are discussed. Some difficulties of the method of zero-range potentials and its possible elimination are discussed. 
We argue that such advanced method of ZRP potential can be applied to realistic electron-molecular processes. 
The method reduces electron-molecule scattering to generalized eigenvalue problem for hermitian matrices and 
admit fast numerical scheme. A noteworthy feature of the method is direct possibility to calculate the wave 
^ ' functions (partial waves). The theory is applied to electron-uracil scattering. Partial phases and cross-sections 

I at low energies are evaluated and plotted. 
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^ ■ A growing interest to electron - big molecules (particularly - DNA and RNA bases) scattering leads to 
a development of one-center expansion of molecular potential methods [1, 2]. However, such expansion 
seems to become non-effective while molecular dimension grows: a molecular potential may have many 
local minima which at the one-center expansion leads to extremely non- smooth one-center radial matrix 
^ potential. Such difficulty does not appear within the method of non-overlapped potentials (NP). 
^ The notion of the non-overlapped potentials (named as the muffin-tin ones) appears in a context of solid 
Q-i, state theory [3] and, further, in scattering theory; it exploits the idea of division of a molecule area into 
^ regions - atomic spheres, which surround each atom, such that originally these spheres touch but do not 
overlap. Inside each of the atomic spheres we replace exact potential by a spherical averaged. Outside 
the spheres, we can replace the potential by a constant potential (ordinary taken for convenience to be 
5^ ' zero). This concept is useful in treatment of both polyatomic molecules and solids. For molecules, this 
approximation gives a model for calculating electronic energy levels [4] and scattering phases [5]. For 
crystals the application of the Bloch conditions instead of the boundary condition of scattering theory leads 
directly to the Korringa-Kohn-Rostoker formalism [6]. 

There is a special case of the non-overlapped potentials in which radius of the atomic spheres is zero - 
zero-range potentials (ZRP). The method of ZRP was proposed in 1936 by E. Fermi [7] and has wide appli- 
cations in photodetachment [8], inelastic scattering [9], and other problem of the quantum physics [10]. In 
such applications a potential is represented by a boundary condition on a wave function and acts to s-states. 
There are also generalizations to higher orbital angular momentum states [11, 12]. To our knowledge none 
of these generalizations have not applied to real molecular processes. It is interesting to characterize some 
of the sources of the method disability: (i) in realistic situation a choice of the potentials parameters is very 
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complicated (ii) eigenfunctions for ZRP are singular and sometimes it is not clear how to avoid ambiguity 
in the treatment of this singularity (iii) it is not clear how to improve results of computations when neces- 
sary. In this paper we show that ZRP method (without adjusted parameters), which include some positive 
features of NP, can cope with realistic processes. 

In the section n we show links between zero-range and non-overlapped potentials and discuss how to 
avoid the problems. We argue that the parameters in generalized boundary condition at an atom can be 
obtained by solving one-dimension Schrodinger equation. If parameters is known, then partial phases and 
cross-sections is readily obtained via matrix eigenvalue problem [12]. 

In the section III we consider an application - low energy electron-uracil (U) scattering. We explain 
how to extract ZRP parameters from the information provided by a quantum chemistry package. We also 
summarize the computational steps involved in our calculation and present results of our calculations: 
partial phases, total and partial cross-sections. 

We conclude with a short summary and perspectives. 

II. THEORY 

As it was mentioned in Introduction, we will consider a polyatomic molecule as a system of spherically 
symmetric NPs, which we denote as Vj(|r — aj|). It is convenient to write the partial wave \l'(r) as the 
linear combination 

N oo I 

^(r) = ^ ^ Xii^<l>iimir - Sii) (1) 

i=l 1=0 m=—l 

of the atomic waves $j;m(r) 

^iim{r) = i^ipiim{r)Yimir), 

where i,l,m- atom number and angular momentum quantum numbers correspondingly; tpum - real func- 
tions of the radial variable r; lim(r) - spherical harmonics [14]; aj - the vector that marks the point of 
expansion for each atom; k - radial wave number; N - number of atoms in the molecule. The coefficient 
Xiim controls the contribution of the wave to the partial wave. We suppose that outside the potential 
Vi{r) the basis functions {$iim(r)} satisfy the Helmholtz equation [15]. Then outside the NP action the 
function ifjum must be the linear combination C2/ip'* (kr) — Cih\^^ (kr) of the outgoing and ingoing spherical 
Bessel functions [16]. For a partial wave, the coefficients ci,2 differ in the phase factor exp(2i5), where S 
is a phase shift. The mentioned choice reflects the condition of equality of the fluxes of the ingoing and 
outgoing waves at infinity. Without loss of generality we may assume that the function i/jum is normalized, 
so that at infinity 

ipiimir) = COS 6 ji{kr) - sin 6 yi{kr), (2) 

where (we used property 10.1.1 of Ref. [16]) ji{x) and yi{x) are regular and singular spherical Bessel 
functions. To derive differential equation for the functions i/Jumir) we expand the partial wave at the 
vicinity of the point a, in the form 

i' i'ipiimir) + Wiimji{kr)) y/m(r), 

Im 

where r/ = r + a^, and substitute the result to Schrodinger equation. The result can be written as inhomo- 
geneous Schrodinger equation for the function ifjum 

{h - \k'^^ tpiim{r) = -WiimVi{r)ji{kr), (3) 
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where the coefficients wum describe the contribution from waves that come from the other NPs, and hamil- 
tonian H is given by the equation 

2 (^dr^ ~^ r dr ^ + «( ) 




FIG. 1: Electrostatic potential for ground state of uracil as function of x,y. Potential was computed in ab initio calculation and plotted in 
Hatree units. 

Using the properties of the Helmholtz equation solutions [14] one can show that the coefficients wum 
satisfy the equation 

><Q\lZirnicos6 - sin 6 yi,{k\^,,\)) , (4) 

where we denote a.ij = ai—aj ; symbol Qll^i^ denotes the integral of three spherical harmonics, multiplied 
by Att. The upper indices mean that the spherical function with indices h, nii should be complex conjugate. 
This integral can be expressed in terms of Klebsh-Gordon symbols, or, alternatively, via the 3jm Wigner 
symbols [14]. In particular, by the Wigner-Eckart theorem the integral is proportional to Clebsh-Gordan 
coefficients [15]: 



Ohm. _ / (2/2 + l)(2/ + l) ,,0 
'^hmilm — \l ^'^ (2/i + 1) ^lOhO^lmhrnz- 



imi 



We can see that these symbols are real. From the coefficients properties [14] it follows, that the nonzero 
symbols correspond to the selection rule mi = m + m2, therefore the sum over m2 may be performed 
explicitly. While account the big / values one can use the selection rule: I + li + I2 - even. 

Note that, the function ^um{r) does not coincide with the radial component of the partial wave for Vi{r), 
because of the nonzero right-hand side of the equation (3). This means that, in general, molecular partial 
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FIG. 2: The scattering phases Su{E) for different atoms and angular momenta Z = 0, 1 as functions of incident electron energy. Solid lines: 
carbon; dashed lines: oxygen; bold lines: nitrogen; dash-dotted lines: hydrogen. 

shifts 5 does not coincide with partial shifts 5ii for a NP. However, partial phases 5ii, coefficients wum, and 
molecular partial shifts 6 are linked by the equality 

Wiim = sin 5 (cot 6ii — cot 6). 

Combining this expression with (4) we see, that the unknown xum satisfy the system of linear algebraic 
equations, which coincide with result of work [12]. In actual calculations values of / are always restricted. If 
the maximum quantum number for all atoms are the same and equal to L, we have the system of N(L + 1)^ 
equations, which gives the same number of partial waves. 

In practice it is convenient to work with matrices and column vectors. Denote the column vector com- 
posed of the coefficients xum as Next define the matrix S, so that 

h+l 

^ilm,jhmi = ^^Qllrn2lmjl2m^ij\)^l2m,2{^ij)^ (5) 

l2 = \h-l\ 

where m2 = mi — m, and matrix N 

h+l 

Klm,jhm^= Yl '''"QhZlrnyhikM)Yi,rn2i^ij) (6) 

l2 = \h~l\ 

for the case when i and j are different, and in the case i = j nonzero elements are given by Num^um = 
cot 5ii. The numbers Urn and jhmi play the role of the multi-indices. It is seen from the last expressions 
that, the matrix S is positive definite and symmetric with respect to transposition of jhmi and Urn, made 
simultaneously with the operation of complex conjugation (hermicity). The matrix N is also hermitian. The 
system for the unknown in the matrix form looks as the eigenvalue problem 

N|x)=AS|a;) (7) 

for the eigenvalues A = cot 5, which are real. 

A result of the eigenvalue problem (7) solution is a set of vectors \x) and the partial phases 5. If the 
molecule has a symmetry, it is possible to classify the phases by irreducible representations of the symmetry 
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TABLE I: In the second column corresponding atomic radii (di) are presented. In the third and fourth columns Xi,yi are coordinates of the 
position vectors a^. We also assume Zi = 0. 
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1.23 
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H7 


0.64 


-0.03767 


-5.22343 


H8 


0.64 


4.08399 


-2.98680 


H9 


0.64 


-3.78218 


-2.82664 


HIO 


0.64 


-0.04204 


3.74037 


Oil 


0.93 


-4.30756 


1.79697 


OI2 


0.93 


4.27587 


2.02817 



group. The multiplicity of degeneration coincide with the dimension of the irreducible subspace. Note, 
that partial phases, assumed as functions of energy, only of different symmetry can intersect, while the 
intersection of the partial phases of like symmetry is impossible. 
Using known vectors |x) we can calculate the partial harmonics 

A{n) = '^Xiimexp{-iknai)Yira{n). 

ilm 

Partial harmonics of the molecular partial waves play the role of spherical harmonics for the spherically 
symmetric potentials. In analogy with the spherical harmonics, partial harmonics are normalized with 
respect to integration over angles to the unit. This means that vector |x) must be normalized, so that 
(x|S|x) = 1. The partial harmonics also may be applied to calculation of the expectation value of an- 
gular momentum operator or angular distribution of scattering for any partial wave. For example, partial 
differential cross section has the form [17] 

The second step of the method is, if necessary, an evaluation of the partial wave for the system of NPs. 
For this purpose the equation (3) is integrated numerically with account the behavior at infinity (2), which 
is already known. After that the partial wave is written as the linear combination (1) of atomic waves 
$i;m(r). Further the perturbation theory can be applied. Methods of perturbation theory (for example MP2 
[18]) were found good for molecular states calculations, it allows to expect to achieve good results by rather 
simple scheme. 

III. EXAMPLE OF URACIL 

As an illustration, we consider low energy electron-molecule scattering for uracil molecules. The com- 
putational steps involved in our calculation can be summarized as follows. 

(i) Choice of target basis and eigenstates. We restrict our consideration to the fixed-nuclear approxi- 
mation and assume that ground molecular state remains unperturbed during the scattering. The single- 
determinant molecular wave function (in basis set 6-31G(d)) and static potential were obtained in ab initio 
calculation with quantum chemistry package. 
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(ii) Choice of effective radii. To compute NPs we divide all space to non-overlapping areas |r — aj| ^ di 
corresponding to separate atoms with effective radii rfj. In order to evaluate the atomic radii d-i we used the 
bond lengths, which were obtained using a geometry optimization procedure, and the following equations 

d^i ~\~ dij ^ I ci.j I * 

The effective radii are listed in Table I. 

(iii) Choice of NPs. To avoid the computational complexity related with coupling terms, we average 
static interaction over angles and compute the radial NPs which were expressed in the following form 

V;t(r) = r ^ d 

r 

where q{r) denote effective charges. Our calculations show that q{r) is smooth function. Calculated 
static potential are plotted in Figure 1 . To include the nonlocal exchange interaction we used local density 
approximation, originally proposed by Slater [19, 20] for atomic problems. Further the approximation 
of the local potential was applied to scattering problems [21, 22, 23] and calculating molecular quantum 
defects for small closed-shell target molecules [24]. In the paper of Hara [23] the following effective 
potential, depending on energy of impact electron is introduced 

VUr) = --Mr) fi + ^ln 

where kj!{r) = [37r^p(r)]^/'^ - Fermi momentum, and p(r) denotes averaged (over angles around the atom) 
electron density. These densities were obtained in ab initio calculation via density matrix. The parameter 
7] is given by the equality 

r/= [P + 2/ + A;2(r)]i/V^F(r), (9) 

where / - molecule ionization energy (in a.u.). The equation (9) emerges from the assumption that the scat- 
tered electron and the electron in the highest energy bound state (i.e. Fermi electron, which has momentum 
fcF(r)) move in the same potential field. For the uracil molecule the experimental value of the ionization 
energy is estimated as / = 8.35 eV. The exchange Hara's potential (8) is attractive (at arbitrary k), because 
the expression in brackets is always positive. When the scattering electron energy (A;^/2) grows, the ex- 
pression at the bracket uniformly decrease, that qualitatively in accordance with the these about decreasing 
with energy contribution of exchange integral. 

(iv) Computation of partial shifts. To obtain the partial shifts 5 a, we solved numerically the radial 
Schrodinger equation with the potentials, which were obtained from step (iii). The partial shifts were 
calculated for angular momenta / = 0, 1 and are ploted in Figure 2. 

(v) Computation of the matrices S and N. These matrixes (for every irreducible representations) were 
constructed by the formulas (5) and (6). 

(vi) Computation of 5 and a. Using the matrix elements assembled from step (v), the molecular partial 
shifts were calculated from eigenvalue problem (2) using Cholesky decomposition. The results of the 
computations at different energies are summarized in Figures 3 and 4. The integral cross section, averaged 
over molecular orientations, were found by formula a = {ATi/k"^) ^ sin^ 5. The integral cross sections are 
plotted in Figure 5. 

(vii) Computation of resonances. The resonance positions and widths were found via the equations 

5{E) = n/2, r = 26'{Ey^ 



1+V 



(8) 



and are listed in Table II. 
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FIG. 3: Computed partial phases 5(E) for A' symmetry as functions of incident electron energy. 
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FIG. 4: Computed partial pliases S{E) for A" symmetry as functions of incident electron energy. 

The Figure 5 shows that resonance positions (2.16, 5.16, 7.8 eV) found in i?-matrix calculations [2] agree 
with our calculations. The sharp peak at 2.2 eV in Figure 5 may also arise from long-lived anion state, 
whose lifetime is sufficiently long to allow nuclear motion and, as consequence, possible fragmentation. 
Indeed, recent experiments [25, 26] showed effective destruction of uracil through dissociative electron 
attachment (DEA). Also a large peak around 1.0 eV was observed. In Ref. [26] authors explain that the 
sharp peak, in their DEA measurements, can be associated with vibrational Feshbach resonance, that is, 
exited vibrational levels of the dipole bound anion states of these compounds. In should be noted, that we 
exclude (for simplicity) dipole interaction. Thus, our sharp peak indicates another possible explanation, 
which, however, does not exclude the arguments of Ref. [26]. 
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TABLE 11: Resonance parameters. 



Symmetry 


-Ercs 


r 




(eV) 


(eV) 




i.Z 


1 1 




2.2 


0.3 




6.7 


1.5 


A" 


1.9 


4.4 




5.6 






6.1 






FIG. 5: Cross sections as functions of incident electron energy. Solid line: our integral cross sections; upper points: partial cross sections with 
the symmetry A'; lower points: partial cross sections with symmetry A". 



IV. CONCLUSION 



We have discussed the "multiple- scattering" method emerged from the zero-range and muffin-tin theo- 
ries. The method reduces electron-molecule scattering to generalized eigenvalue problem (7) for hermitian 
matrices. A noteworthy feature of the method is direct possibility to calculate the wave functions (par- 
tial waves). This way is based on one-dimensional inhomogeneous Schrodinger equation (3) for radial 
component ipumir) and only requires the atomic potential and partial phase of given partial wave. Thus, 
this approach may be especially important in relation to Bardsley-Fano theory [27, 28] and dissociative 
attachment at low energies. 

The preliminary results on uracil presented here are quite promising. Particularly, our positions of the 
shape resonances demonstrate good agreement with the i?-matrix calculations [2]. Thus, this method can 
be applied to another DNA bases, such as adenine, cytosine, thymine etc. In this connection, it is interesting 
and important to extend the method to the case of long range interactions, such as dipole interactions. Note, 
this method gives fast numerical scheme and, therefore, we have good base for future approximations and 
modeling of electron-DNA scattering in sense of the works [29]. 
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